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Introduction 

Turbulent flames are often envisioned as an ensemble of random vortices interacting with the combustion 
process. A better understanding of the vortex-flame interactions therefore would be useful in improving the 
modeling of turbulent diffusion flames. Substantial simplification may be made by investigating controlled 
interactions in a laminar flame, as opposed to random interactions in a turbulent flame [1]. The general goals of 
the research project are to improve our understanding of (A) the influence of buoyancy on co-flow diffusion 
flames, and (B) the effects of buoyancy on vortex-flame interactions in co-flow diffusion flames. As a first step 
toward Objective (B), we conducted a joint experimental and numerical investigation of the vortex-flame 
interaction. Vortices were produced by mechanically pulsing the fuel flow at a low frequency, e.g., 10 Hz. 

Experiments were conducted using a non-flickering Burke-Schumann [2] flame in both microgravity (pg) and 
normal gravity (lg) as a means of varying the buoyant force without modification of the pressure (i.e., density). 
The effects of buoyant convection may then be determined by a comparison of the pg and lg results. The pg 
results may also reveal the important mechanisms which are masked or overwhelmed by buoyant convection in 
lg. A numerical investigation was conducted using a validated, time-accurate numerical code to study the 
underlying physics during the flame interaction and to assist the interpretation of the experimental results. 

Numerical Model 

The mathematical formulation was based on the Burke-Schumann flame sheet model [2]. The time-dependent, 
axisymmetric Navier-Stokes equations coupled with the continuity and mixture fraction equations were solved 
by a semi-implicit scheme, which was modified from that reported in [3], by application of a flux-corrected 
transport (FCT) term [4] to the QUICK-based convective flux and a second-order central differencing to the 
diffusional flux. The projection method was retained to solve the velocity and pressure coupling through a 
direct solver. Two computational grid distributions, 180 (axial) by 80 (radial) and 120 by 60, were employed to 
cover a computational domain of 20.1 cm (axial) by 2.33 cm (radial), matching the chimney of the experimental 
burner. The fuel jet diameter was set to 0.56 cm, matching the inner diameter of the experimental burner. The 
two different grid distributions predicted the same dynamical characteristics but the denser grid produced a 
better resolution for isotherm and vorticity images. Steady-state solutions were obtained for the non-pulsed 
case and used as the initial condition for the simulation. The computation was done using a computer 
Workstation (RISC-based CPU-50 MHz) or a Personal Computer (486-66 MHz, DX2). With the 120 by 60 grids, 
the CPU times for each time step (0. 2 ms) are 5.1 s and 10.5 s for the Workstation and PC, respectively. 
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Modeling Results 


The BSDF studied in this paper does not exhibit flickering in lg (normal gravity) or Og (zero gravity) for both 
the simulation and experimental results. In Og, the predicted flame (stoichiometric boundary) is wider and 
slightly shorter than in lg , in qualitative agreement with [5]; the luminous region of experimental pg diffusion 
flames were slightly taller than the lg flames, however. To simulate the experimental pulsation of the fuel jet 
due to opening and closing of a solenoid valve, a sine-wave forcing with an amplitude of 80% of the mean fuel- 
jet velocity (up was imposed on the fuel jet inlet boundary condition. The simulation shows that the imposed 
sine-wave perturbation results in periodic flame oscillation, responding to the imposed forcing and with a 
phase lag of between the maximum flame height and the maximum Uj during one forcing cycle. No tip cutting is 
seen at low velocity, e.g., uj = u a (coflow air velocity) = 5 cm/s. The cutting occurs at higher velocities. 

When the fuel-jet velocity is increased, flame-tip cutting (i.e., pinching and separation of the flame tip) is 
predicted for uj = u a = 18 cm/s, in both Og and lg, cf. Fig. 1. The flame contours at ten phase angles separated by 
10 ms (or phase angle of 36 degrees) are shown. The 18 cm/s emulates the non-perturbed 20 cm/s condition in test 
3 (cf.. Table I). The arbitrarily chosen reference time, i.e., t = 0 ms, has a nearly matched flame shape at t = 50 
ms for the Og and lg flames. The flame-tip cutting for the lg flame occurs at 10 ms < t < 20 ms, whereas in Og, the 
flame tip detaches at 20 ms < t < 30 ms. Unlike the Og flame, the detached flame parcel in lg disappeared 
within the 10 ms interval. The flame-tip cutting is a result of an enhanced local mixing which leads to a local 
mixture fraction below the stoichiometric value. This enhanced mixing can be attributed to the vortex-flame 
interaction in a lg pulsed methane-air diffusion flame [6], in which the vortex was formed due to the buoyancy 
acceleration. The present simulation and experiments showed a similar tip-cutting in BSDF in jxg. The 
predicted flow-field isotherm, vorticity, and velocity contours revealed a large vortex outside the flame in the 
BSDF with Uj = u a = 18 cm/ s for both the Og and lg conditions. The pulsation of the fuel jet was responsible for 
the formation of an initial vortex near the fuel nozzle. The tip cutting was seen to coincide with the formation 
of a new vortex near the flame tip in lg, and below the instantaneous half flame-height in 0g. 

Experiment 

Tests in the 2.2-second Drop Tower at the NASA Lewis Research Center [ 7 ] were conducted by pulsing the fuel 
flow in a laminar BSDF. The nitrogen-diluted fuel (CH4/ N2 50% by mole) was used to reduce the sooting since 
sooting obscures the true flame shape as determined from the video image. Dry air was the oxidizer, and all 
tests were performed at approximately 0.98 atm. The small cylindrical burner consisted of a glass chimney (i.d. 
of 4.66 cm and length of 20 cm) surrounding a thin-walled (0.04 cm) stainless-steel tube (i.d. of 0.56 cm), with a 
fuel-to-oxidizer tube diameter ratio of 0.12 [3]. Flow restrictions and straighteners were used to produce a flat 
velocity profile at the entrance of the annulus. The fuel tube had a straight length with a length-to-i.d. ratio 
of approximately 27, presumably having a parabolic flow at the exit. A single layer of fine wire mesh was 
placed across the outlet of the chimney to keep sand (the drop tower deceleration material) out of the burner. 

The tests were performed with uj in the range of 10 to 30 cm/s, measured with a mass flowmeter and without the 
pulsation. The fuel flow was pulsed at a low frequency in the range of 6.7 to 12.5 Hz by opening and closing a 
solenoid valve (specified response time of 5-10 ms) in the fuel supply line. The solenoid valve was closed for 
approximately 10 ms and then held open for the remainder of the cycle. Two particle filters (valve coefficient 
of 0.15) were positioned between the solenoid valve and the burner inlet to produce a small pressure drop, and 
thus limit the amplitude of the pulsating inlet velocity. The annular air velocity was nominally set to be equal 
to the non-perturbed fuel velocity. Regrettably, after cessation of these tests (due to a shutdown of the Drop 
Tower for rehabilitation) a noticeable leak in the oxidizer system between the flow meter and the burner was 
discovered. Even with the leak, the air velocity was such that the flame did not flicker when the fuel flow 
was not perturbed. 

The flames were ignited in lg prior to the drop, with a hot-wire ignitor which was retracted (about 1 cm) out of 
the flame following ignition. In a typical test, the flame was allowed to develop without pulsations for 2 s, 
and then pulsed for 3 s in lg. After a pulsation-free period of 5 to 7 s, the flame was dropped and the pulsations 
were resumed (at 0.1 s into the drop) for 2 s in pg. Although lg ignition leads to residual buoyant flows, the gas 
velocities were such that buoyant effects should be convected away in less than 1 s (e.g., a flame-zone velocity 
of 40 cm/s would convect buoyant effects from the 20-cm chimney in 0.5 s). Furthermore, the experiments of [8] 
showed that flames ignited in lg will more quickly reach (nearly) steady conditions during the drop test, than 
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flames ignited in pg. 


For the present investigation, the instrumentation was strictly photographic with a color video camera 
recording the flame's side view at 30 frames per second (or 60 fields per second) [3]. The flame response to the 
flow pulsations was quantified by a measurement of the flame height as a function of time, using an image- 
processing and object-tracking workstation at NASA Lewis [9]. When the tip cutting occurred, the detached 
flame parcel was included in the flame height determination [3]. 


Table I: Summary of Test Results. 



Test Conditions 

Flame Response, <(>, 

Flame Height (cm) 



Forcing 

Fuel 

(and tip cutting) 


Mean (and std. dev.) 


Frequency 

Velocity 



without 

with pulsations 

Test 

(Hz) 

(cm/s) 

lg 

MK 

pulsations 


UK .... _.. 

1 

10 


l(n) 

1 (n) 

1.2 

0.9 (0.1) 

0.7 (0.1) 

2 

10 

15 

3 (n) 

Kn) 

2.1 

1.9 (1.1) 

E — 

3 

10 

20 

2(y) 

2(n) 
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4 

10 

25 

2(y) 

2 (y) 

3.4 

4.1 (2.6) 

3.3 (1.3) 

5 

10 


i (y) 

2 (y) 

4.3 

5.0 (2.1) 

4.2 (1.8) 

6 

6.7 

20 

2 (n) 

1 (n) 

2.6 

2.4 (1.5) 


7 

8.3 

20 

2 (n) 

1 (n) 

2.6 

2.6 (1.7) 

2.3 (0.6) 

3 

10.0 

20 

2 (y) 

2 (n) 

2.7 

2.9 (2.0) 

2.3 (0.7) 

8 

12.5 

20 

2 (y) 

3 (y) 

2.7 

2.5 (1.8) 

2.1 (1.0) 


Experimental Results 

The video recording showed that the names puisatea in response to the flow perturbations, and tip cutting was 
sometimes observed. When the flow was pulsed, the flames appeared to be strongly flickering, although the 
oscillations were the result of the forced perturbations and not buoyancy-generated instabilities. In addition to 
the fundamental frequency (i.e., showing the perturbation frequency in the response), the flame also exhibited 
a subharmonic response which varied with the gravity level, frequency, and flow rate. The subharmonic flame 
response can be seen in the photographs reported in [3]. For example, when the fuel jet with a mean velocity of 
18 cm/s and pulsed at 10 Hz (as in Fig. 1), the flame responded subharmonically at 5 Hz in both pg and lg as 
shown by the flame-height plots in Fig. 2; the flame-tip cutting was observed in lg but not in pg. In other tests, 
however, the flame responded differently in pg and in lg. A summary of the test results is shown in Table I. 

A flame frequency response parameter, <j>, can be defined to be the ratio of (number of fuel-flow-cycles imposed) 
to (number of flame-shape-oscillations in response ) to describe the dynamic characteristics. For example, <|> is 
equal to 2 for both gravity levels shown in Fig. 2. As seen in Table I, the mean flame heights were slightly 
greater in lg than pg, in contrary to the steady-state predictions of [10-12] or earlier observations of non- 
perturbed BSDF [5, 8, 10, 11] in which the pg flames were slightly taller than the lg flames. However, the 
observation of the mean flame height being slightly less in pg is in agreement with the numerical predictions 
presented in this paper. 

The flame height as a function of fuel flow is shown in Fig. 3. For both lg and pg conditions, the mean flame 
height increased linearly proportional to the fuel flow rate (as measured when it was not perturbed). The 
perturbed flows were not measured; thus, the figures can be misleading because the mean perturbed flow is lower 
than the non-perturbed value. For this reason, the height of the non-perturbed lg flame was included in the 
figures. In general, the flame-height oscillation amplitude increased with increasing flow rate, and the 
relative position of the mean height within the amplitude range appeared to be related to <|>. 

The variation in forcing frequency over the range of 6.7 to 12.5 Hz had a limited effect on the flame response in 
lg. All four flames in Table I (i.e. with four different uj values) responded bimodally (i.e., <(>ig = 2), and since 
the fuel flow was approximately fixed, the flame height oscillations were also fairly similar in the observed 
maximum, mean, and minimum heights. The test procedure resulted in a reduction in fuel flow with increased 
frequency, because the flow rate was set under no perturbations and a constant valve closure time for each cycle 
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was used for all frequencies (i.e., 10 ms). The most significant effect of the forcing frequency was on the tip 
cutting, which was apparent in the video record for the high frequencies but not the low frequencies. However, 
it was also possible that the flame tip cut between video frames because the tip often pinched partially when 
cutting was not apparent. In pg, the flame response becomes increasingly subharmonic as the forcing frequency 
was increased. Other than the change in <|>, the height oscillation of the pg flames was only moderately 
affected by the forcing frequency (with test 8 at <tw = 3 showing the most dramatic change, cf. Table I). The 
flame frequency response, <|>, generally increased with increasing fuel flow and forcing frequency in pg; but <f> 
varied non-monotonically with fuel flow and it was not affected by the forcing frequency in lg. However, the 
tip cutting occurrence increased in both lg and pg with increasing fuel flow or forcing frequency. 

The predicted flame heights at different fuel-jet velocities are also shown in Fig. 3. The predicted mean flame 
height increased linearly with increasing us, similar to the trend observed in the experiments. The predicted 
mean flame heights were longer than the experimental data, but the discrepancy decreased as the fuel-jet 
velocity was increased. The predicted mean flame length regresses to zero at u; = 0 while there is an offset for 
the experimental data. The offset is due to the perturbation technique which diminished the flow rate from 
that reported. The predictions also showed <j> = 1 for all the conditions examined; the experiments, however, 
had different <> values, cf. Table I. 

Discussion 

The flow characteristics during the tip cutting can be deduced from numerical visualization based on the 
predicted isotherms and vorticity. For example, in Og and lg flames, large vortical structures (or concentrated 
vorticities) were formed outside the stoichiometric surface. These concentrated vorticities were caused by the 
sharp velocity gradient adjacent to the injector wall, and first appeared near the fuel nozzle and were then 
convected downstream. In Og, the vorticity strength decayed rapidly when the vortex was convected 
downstream, but this rapid decay was not observed in lg. In lg, the downstream maximum-vorticity was found 
to be comparable to (or even slightly higher than) the maximum at the fuel jet exit. Both the buoyancy 
acceleration and volumetric expansion significantly affected the species transport (i.e., the mixture fraction). 
For e xam ple, the lg flame length continued to increase beyond the phase angle or the time instant of the 
maximum fuel-jet velocity at nozzle exit because of buoyant convection of the stoichiometric surface. The 
volumetric expansion, however, appeared to "hold up" the flame tip at certain phase angles during the forcing. 
In Og, the predicted tip cutting was located at below half of the instantaneous flame height at a phase angle 
having a zero fuel-jet velocity at the nozzle exit. In lg, the tip cutting occurred near the flame tip at a phase 
angle before the zero fuel-jet velocity was reached. The unsteady flow due to the pulsation of the fuel jet was 
responsible for the high vorticity value at the nozzle exit and responsible for the vortex-flame interaction. 

The buoyant acceleration in lg enhances the downstream convection of combustion products, and the vortex 
growth. In Og, the volumetric expansion in the fuel-lean region appeared to be a sink term in the vorticity 
budget; whereas in lg, the buoyant acceleration was a source term in the fuel-lean region. The observation was 
consistent with theoretical analysis of the vorticity equation of a jet diffusion flame reported in [13]. 


The numerical simulation predicted the BSDF tip-cutting in response to the periodic forcing of the fuel jet. The 
simulation, however, did not predict the experimentally observed bimodal or trimodal response. On the other 
hand, the current numerical model was able to predict the subharmonic frequency (resulting from vortex 
merging) of the flame flicker for the condition reported in [14]. Thus, a parametric study was conducted to assess 
such effects on the flame response as different velocities, fuel-jet velocity (spatial) profile, driving waveform 
and driving frequency. Decreasing u a to half or one-third of the fuel-jet velocity (i.e., to 9 or 6 cm/ s) was used to 
simulate the leak in the experiments, but it yields only a slight decrease in flame length and the tip cutting 
was still predicted. The predicted $ value, however, did not change. Using a uniform or parabolic velocity 
profile at the fuel-nozzle exit did not change the prediction of <(> either. A slightly longer maximum flame- 
height was predicted when the parabolic profile was used. 

Three waveforms were used as the fuel-jet boundary condition in an effort to evaluate the effects due to the 
perturbed fuel-jet velocity profile since the profile was not quantified in the experiments. None of these 
waveforms, however, predicted a bimodal or trimodal response, nor a 4> * 1 response. The waveforms evaluated 
included (A) using a zero velocity for the initial 10 ms and followed by 20 cm/s for the remaining 90 ms of a 100 
ms forcing cycle, (B) a zero velocity for the initial 10 ms, 32.5 cm/s for next 5 ms, and 19.27 cm/ s for the 
remaining 85 ms, and (C) a zero velocity for the initial 10 ms, 53 cm/s at 10 ms and followed by an exponential 
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decay to preserve the averaged flow rate. The predicted flame length varied with these different waveforms; 
for example, the tip-cutting was not predicted. It was also observed that higher harmonic frequencies in the 
flow field were predicted with a 1-Hz or 2-Hz forcing of the fuel jet, when a higher-order scheme was used 
(i.e., the current FCT scheme). This higher harmonic frequency response, however, was not predicted when a 
lower-order scheme (e.g., an upwind) when the same grid resolution was used. Thus, a higher-order scheme is 
required for accurate prediction of dynamic response of the chemically reactive flow investigated in this paper. 

To assess the effects due to the model assumptions on the prediction of flame response (or the failure in the 
prediction of the bimodal or trimodal response), the finite-rate chemistry, equal diffusion, and unity Lewis 
number assumptions were relaxed. The simulation based on the one-step global kinetics of Westbrook and Dryer 
[15], differential diffusion of different species, and with the non-unity Lewis number effect accounted for in the 
energy equation predicted essentially the same flame response, i.e., only with $ = 1. The results suggested that 
the flame dynamics were dominated by convective flow, and the transport properties and finite rate chemistry 
played a secondary role. It should also be mentioned that the differential pressure across the solenoid valve 
was measured and there was no evidence to suggest that a sub-harmonic response existed in the valve opening 
and closing . Therefore, the failure to predict the bimodal or trimodal response was not due to the numerical 
model assumptions, nor was it due to the presence of a subharmonic response in valve characteristics. Further 
investigation is needed to resolve the discrepancy. 
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Figure 2. Observed flame height vs. time in lg and jig with a 10-Hz driving (test 3, Table I). 
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Figure 3. Flame height vs. non-perturbed fuel-jet inlet velocity; detached flame parcels are included and <)> 
values are shown for fuel-jet inlet velocities of 5, 10, 20, 25 and 30 cm/ s. 
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